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Abstract 

Within a perturbative approach we investigate decays of charmonium states 
into baryon-antibaryon pairs. Using a recently proposed wave function for 
the nucleon and suitable generalizations of it to the hyperons and decouplet 
baryons, we obtain the decay widths for the BB channels in reasonable 
agreement with data. An important difference to previous work is the use of 
the c-quark mass in the perturbative calculation instead of the charmonium 
mass. As a consequence of this feature our approach possesses, the J /if) 
and the if>' decay widths do not scale with a high power of the ratio of their 
masses. 
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1 Introduction 



The decay J ftp —>■ pp has been investigated within perturbative QCD by Brodsky and 
Lepage first Later on this analysis has been repeated several times, e.g. 0, [J, and 
even extended to the AA decay channel. It has been argued that the dominant dynamical 
mechanism is cc annihilation into three gluons and subsequent creation of light quark- 
antiquark pairs forming in turn the final state baryons. Three is the minimal number of 
gluons allowed in J/ip decays; cc annihilations through one or two gluons are forbidden 
by colour and C-parity, respectively. Contributions from annihilations through more than 
three gluons constitute higher order corrections. The dominance of annihilation through 
gluons is most strikingly reflected in the narrow width for hadronic channels in a mass 
region where strong decays have typically widths of hundreds of MeV Q . The dynamical 
suppression at work here is perturbative QCD of higher orders (the total hadronic J/tp 
decay, for instance, is an O(af) process) and is customarily regarded as evidence for the 
Zweig rule. Since the c and the c quarks only annihilate if they are separated by distances 
less than about l/m c and since the average virtuality of the gluons is about 1 GeV 2 
one may expect perturbative QCD to be at work although corrections are presumably 
substantial. Indeed as the previous perturbative analyses, performed in the standard 
hard scattering approach, i.e. in collinear approximation, showed the J/ip decay into pp 
seems to be fairly well described. This is in marked contrast to the case of the nucleon 
form factor where soft physics seems to dominate in the experimentally accessible region 
of momentum transfer [||, |7J. A point to criticize in these studies of the J/ip decays 
is the treatment of the strong coupling constant a s . Since, as we mentioned, the average 
virtuality of the gluons is about 1 GeV 2 one would expect a s to be of the order of 0.4 
to 0.5 rather than 0.2 to 0.3 as is customarily chosen p, Since a s enters to the sixth 
power into the expression for the width a variation of a s from, say, 0.2 to 0.3 would lead 
to a change by a factor of 11 for the width. Thus, a large factor of uncertainty is hidden 
in these calculations preventing any severe test of the wave function utilized. 

In constrast to previous works |], |^, [3| we will not use the collinear approximation but 
rather the modified perturbative approach of Sterman et al. || in which transverse degrees 
of freedom are retained and Sudakov suppressions, comprising those gluonic radiative 
corrections not included in the evolution of the wave function, are taken into account. An 
important advantage of the modified perturbative approach is that the strong coupling 
constant can be used with a renormalization scale depending on the momentum fractions 
the quarks carry und thus large logs from higher orders of perturbative QCD are avoided. 
This choice of the renormalizaton scale entails singularities of a s which are, however, 
compensated by the Sudakov factor. Hence, there is no uncertainty in the use of a s . This 
is to be contrasted with the standard perturbative approach where either a s is evaluated 
at a renormalization scale that is a constant fraction of M 2 , or at a momentum fraction 
dependent scale in which case a s is to be "frozen" at a certain value (typically 0.5) in 
order to avoid uncompensated a s singularities in the end-point regions. The modified 
perturbative approach possesses another interesting feature: the soft end-point regions 
are strongly suppressed. Therefore, the bulk of the perturbative contribution comes from 
regions where the internal quarks and gluons are far off-shell. In contrast to the nucleon 
form factor the J /if) — > BB amplitude is not end-point sensitive. The suppression of the 
end-point regions does not, therefore, lead to a substantial reduction of the J/tj) — > BB 
amplitude. For the same reason, the size of that amplitude does not exhibit an extreme 
sensitivity to the baryon wave function utilized in the calculation as is, for instance, the 



case for the baryon form factor. 

In the calculation of the decay widths we will make use of a (valence Fock state) 
wave function for the nucleon that we proposed recently ||. That wave function was 
constructed in the following way: In accord with the findings reported in |J where it 
was shown that a reliably calculated perturbative contribution to the nucleon form factor 
is very small, the nucleon wave function is demanded to describe the form factor via 
the Drell-Yan overlap contribution for momentum transfers around 10 GeV 2 and to be 
compatible with the available valence quark distribution functions of the nucleon. As a 
third constraint on the nucleon wave function the decay J/ip — » pp was employed. The 
nucleon wave function proposed in || will be suitably generalized to the case of hyperons 
and decuplet baryons. The modified perturbative approach is then used to calculate the 
widths for the J/ip decays into pairs of octet, B 8 , and decuplet, B w , baryons. We are 
also going to calculate the baryonic decays of the ip' for which recently the first, still 
preliminary, experimental results for channels other than pp have been reported ||. 

The paper is organized as follows: In section 2 we introduce the octet baryon wave 
functions used in the analysis described in detail subsequently and briefly recapitulate a 
few properties of light-cone wave functions. In section 3 we present the wave functions 
for the decuplet baryons. Section 4 is devoted to the calculation of the J /if) — > B 8 B 8 
decay widths. This analysis is extended to the decays into decuplet baryons (section 5) 
and to decays of other quarkonia into BB pairs (section 6). Finally, section 7 contains 
our conclusions. 



2 The wave functions of the octet baryons 

Generalizing the ansatz made for the nucleon in [||, || 10], we write the valence Fock states 



of the lowest lying octet baryons as (the plane waves are omitted for convenience) 

\b 8 ,+) = ^jp / [da;][d 2 kj 1 nxnmi > + 1 nmmx > 



*f 3 8 2 + *SiJ \m.f&f&)) (2-la) 
A , + ) = / [dx] [d 2 kj | < 3 ) - I u a 2 < 2 < 3 ) 

- *23i) K 1 (2-lb) 



where (2.1a) holds for all octet baryons except the A and the E°. Obviously, for the 
proton and the S + fi represents an u quark and f% either a d or a s quark, respectively. 
For the H~ j\ represents a s quark and f2 a d one. The states of the neutron, E~ and 
S° are obtained from those of the proton, S + and H~ states, by exchanging u «-> d, 
respectively. The baryon is assumed to be moving rapidly in the 3-direction. Hence the 
ratio of transverse, kj_j, to longitudinal momenta, Xjp, of the quarks is small and one may 
still use a spinor basis on the light cone. The integration measures are defined by 

[dx] = f[ ^ S(l - E *i) [ d2k J = TTTtW n d2k ^ 5(2) (E k u) • (2.2) 

i=l i l lt)7r ) i=l i 



The quark fa is characterized by the momentum fraction Xi, by the transverse momentum 
h± i as well as by its helicity Aj and colour a». A three-quark state is then given by 

1 



UlAiJ2A 2 J3A3/ 



y/XiX 2 X 3 

The single-quark states are normalized as 



/i 1 ; xi, k ±1 , Ai) | / 2 a2 ; x 2 , k ±2 , A 2 ) | / 3 a3 ; x 3 , k ±3 , A3). 



(2.3) 



4 k x J, X I /f* ; x if k ±4 , A 4 > = 2x J (2 7 r) 3 ^ / ^yA l ^ - a*)^^ - k ±l ). (2.4) 

Since the 3-component of the orbital angular momentum, L 3 , is assumed to be zero the 
quark helicities sum up to the baryon's helicity. (2.1a) is the most general ansatz for the 
L 3 = projection of the three-quark nucleon wave function [IT]. From the permutation 
symmetry between the two u quarks and from the requirement that the three quarks have 
to be coupled in an isospin 1/2 state it follows that there is only one independent scalar 
wave function. If the L 3 ^0 projections are included the entire nucleon state is described 
by three independent functions [11]. In general there are more than one scalar wave 



function for the other octet baryons if SU(3)F-symmetry breaking is taken into account. 
As already expressed in (2.1) we nevertheless assume that each octet baryon is described 
by a single scalar wave function which, for convenience, we write as 

1 TTT^t f^){^p)^ih{ x i ^f) ^(8)0^) kjj • (2.5) 



kj 



We assume that SU(3)p symmetry is only broken by a quark mass dependence of <j) Bs (see 
below). /( 8 ), being related to the wave function at the origin of the configuration space, 
is identified with the nucleon parameter fa whose value was determined in || to amount 
to 6.64 ■ 10~ 3 GeV 2 at the scale of reference /xq = 1 GeV. 

The transverse momentum dependence of the baryon wave function is parameterized 
by a simple symmetric Gaussian 



2\2 
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(2.6) 



and the transverse size parameter arg\ is assumed to be the same for all octet baryons. A 
value of 0.75 GeV -1 is used for that parameter (see 0). 

The last item to be specified in ( |2.5|) is the distribution amplitude, 0^|(x,/ip) = 
4> Bs (xi, Xj, Xk, Hf)i °f an octet baryon B s , which is conventionally normalized to unity 



[dx]0f 2 8 3 (x,/i F ) = 1 



(2.7) 



The distribution amplitude representing the wave function integrated over transverse 
momenta up to the factorization scale, fiF, can be expanded upon the eigenfunctions of 



the evolution kernel being linear combinations of Appell polynomials (see [0, |13|| ) 



n=l 



(2.8) 



where <Pas{ x ) = I2OX1X2X3 is the asymptotic distribution amplitude fl2 |. Evolution is 
incorporated by the factorization scale dependences of /(§) and the expansion coefficients 
B n ■ 



Vln(/i F /A Q c D ), 



s ln(/i F /A Q c D ) J 



(2.9) 
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1 


X\ — X 3 


20/9 


2 


-2 + 3(xi + x 3 ) 


8/3 


3 


2 - 7(xi + x 3 ) + 8(x? + x\) + 4xix 3 


32/9 


4 


(xi - x 3 )(l - 4/3(xi + x 3 ) 


40/9 


5 


2 - 7(xi + x 3 ) + U/3(xl + x\) + 14xix 3 


14/3 



Table 1: Eigenfunctions and reduced anomalous dimensions for helicity 1/2 baryons. 



where /3q = 11— 2/3 n/. The exponents 7™ are the reduced anomalous dimensions. Because 
they are positive fractional numbers increasing with n [12], higher order terms in (|2.8|) 



higher order terms in ( 
are gradually suppressed. The reduced anomalous dimensions and the eigenfunctions 0™ 2 3 
are listed in Tab. [1] where the notation of [Q is adopted. 

In |J the nucleon distribution amplitude was found to have the simple form 



1 + ^23 0*0 + 4^23^) 



60xiX2X 3 [1 + 3xi] 



(2.10) 



The nucleon wave function is fully specified now and, before turning to the discussion 
of the hyperon distribution amplitudes, we note in passing that the probability of the 
nucleon's valence Fock state is 0.17. 

A suitable hyperon distribution amplitude is constructed by taking ( |2.10p and, in order 
to incorporate the empirically known breaking of SU(3)p symmetry, multiplying it with 
a factor 



cxp 



a (8) m s 
Xj 



[2.11) 



whenever the quark j is a strange one. That factor bears resemblance to the BHL expo- 
nential |15| . 



ci(8) is the transverse size parameter already introduced in ( |2.6|) and m s is a 
still to be adjusted parameter related to the strange quark mass. Explicitly our hyperon 
distribution amplitudes read 



0i23(^o) = N S (j)? 23 (x) exp 



°X8) m s 
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N I 

123 1 



x) + 2 



lN 
; 321 



X 



exp 



a (8) m s 



X 3 



feO^o) = N B (f)? 23 (x) exp 
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(2.12a) 

(2.12b) 
(2.12c) 



The constants Nb 8 ensure the correct normalizations (see (|2.7|) ) of the hyperon distribution 
amplitudes (N Bg = 1 for m s = 0). The particular combination of 0at's appearing in the 
A case (2.12b) is required by SU(3)f symmetry. 

In order to take into account evolution properly we expand the distribution amplitudes 
(2. 12a)- (2. 12c) upon the eigenfunctions of the evolution kernel (see (]2.8j) ) up to terms of 
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Set 2 (m s = 150 MeV) 


Set 4 (m s = 480 MeV) 




B 1 


B 2 


B 3 


#4 


B 5 


Bt 


B 2 


B 3 


B, 


B 5 


S 


0.623 


0.284 


-0.085 


-0.285 


0.065 


-0.118 


0.484 


-0.404 


-1.173 


0.358 


A 


-0.360 


0.282 


-0.048 


-0.195 


0.027 


-1.022 


0.478 


-0.182 


-0.650 


0.127 




0.831 


0.201 


-0.083 


0.389 


-0.129 


1.338 


-0.068 


-0.384 


2.943 


-0.775 



Table 2: Expansion coefficients B n (n ) of octet baryon distribution amplitudes for various 
values of the parameter m s . 



order n = 5. In Tab. Q we quote four sets of expansion coefficients B n corresponding to 
the following scenarios: Set 1 is obtained from m s = 0, set 2 from m s = 150 MeV (current 
strange quark mass) and set 4 from m s = 480 MeV (constituent strange quark mass). 
The intermediate set 3 corresponds to m s = 350 MeV ~ (480 2 — 330 2 ) 1 / 2 MeV (difference 
between the squares of strange and light constituent quark masses). This value appears if 
( f2. 11| ) is interpreted as the ratio of the BHL exponentials for a hyperon and the nucleon 
distribution amplitude. 

In Fig. [I] we show contour plots of the four octet baryon distribution amplitudes for 
m s = 480 MeV in order to illustrate the effect of the mass exponential ( |2.11| ). It can be 
seen that, compared to the nucleon case, the maxima of the £ and A distribution ampli- 
tudes are shifted to the right, i.e. to larger x 3 values, whereas that of the S distribution 
amplitude is shifted to the left. On the average, s quarks carry larger momentum fractions 
than d quarks if m s > 0. 

Distribution amplitudes for the octet baryons are also discussed in As compared 
to our ones these QCD sum rule based distribution amplitudes are strongly concentrated 
in the end-point regions and exhibit a more significant breaking of SU(3)p symmetry. 
Moreover, they reveal three pronounced maxima near the end-points. In summary, the 
distribution amplitudes proposed in [16| are very different from ours and also from the 
asymptotic distribution amplitude at experimentally accessible scales. As our distribution 
amplitudes but to a much greater extend, they possess the property that, on the average, 
a u (s) quark in the proton or the S + (H~, A) with the same helicity as its parent baryon 
carries the largest fraction of the baryon momentum. 
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Figure 1: Contour plots of octet baryon distribution amplitudes for m s 
at a scale of 1 GeV. 



480 MeV and 



3 The wave functions of the decuplet baryons 

In analogy to (2.1) we write the valence Fock states of the decuplet baryons A and E* 
with helicity \b w = 1/2 as[] 



|A ++ , + ) 
\B 10 , + ) 



~7T 



|H[d 2 kj | u? ) 



(3.1a) 



f&ffif£)\. (3-lb) 



For the A + (E*~) f\ represents an u (d) quark and / 2 a d (s) one. Again, the states of 
the A , A - and E* + are obtained from those of the A + , A ++ and E*~ by exchanging 
u «-> d, respectively. The generalization to the other helicity states is trivial. There is 
obviously only one independent scalar wave function in the case of the A and we assume 
that the same shall apply to the E*. The scalar wave functions of the decuplet 

baryons Bio (= A, E*) are parameterized in a fashion similar to the octet baryon case: 



24 V2 



&ii£(x,VF) Q(io)(z,kj 



(3.2) 



/( 10 ) is assumed to be equal for all members of the baryon decuplet. We again adopt the 
form (|2.6|) for the transverse momentum dependent part of the wave function, f2(io), with 
the transverse size parameter a^g) replaced by ct(io)- 

1 We omit the discussion of the decuplet states 5(1530) and f2~ since the perturbative approach is 
not applicable to the decays J/ip — > S(1530)S(1530) due to almost zero momentum transfer. The decay 
J /if) — > f2~f2 + is even kinematically forbidden. The decay into X* is the borderline case for the application 
of the perturbative approach. 



Decuplet baryons in helicity 1/2 states have the same eigenfunctions of the evolution 
kernel and the same anomalous dimensions as the octet baryons (see Tab. H). Evolution 
of the helicity 3/2 baryons is, on the other hand, different |T2"|| . We refrain from giving 
details on that case here since we do not consider such baryons. 

In the particular case of the A it is tempting to use a completely permutation symmet- 
ric distribution amplitude since the A is composed of are three light quarks in symmetric 
spin and flavour states. Noting that the permutation symmetric part of the nucleon dis- 
tribution amplitude (|2.10|) is just equal to the asymptotic distribution amplitude, we take 
<Pi2 3 (x) = 4>as(%) and construct s * (x) analogously to s (x) 

= Nx* 0as(x) exp f-^p) • (3-3) 

The evolution of the S* distribution amplitude is treated as in the case of the octet 
baryons by expanding ( |3.3|) upon the eigenfunctions <Pi23{x) up to n = 5. In distinction 
from the octet baryon case, the QCD sum rule based distribution amplitudes for the 
decuplet baryons |17]] are not very different from the ones we are proposing. 



4 Decays of the J/ip into octet baryons 

Now, with the model wave functions at hand, we can calculate the decay width of J/ip into 
octet baryon-antibaryon pairs within the modified perturbative approach. The helicity 
amplitudes of these processes may be decomposed covariantly to read 

Wp 2 ,A 2 K(A) (4.1) 

where ms 8 is the mass of an octet baryon. pi (p?) and Ai (A2) denote the momentum and 
the helicity of an octet baryon (antibaryon), respectively. ub 8 and vb s are their spinors 
(normalized as ub s ub s = 2mg 8 ) and e is the polarization vector of the J/ip- in a leading 
twist perturbative approach the helicity amplitudes are determined by the invariant B Ba 
(C Bs = 0). Implicitly this ensures hadronic helicity conservation. The J/ip — > B$B S decay 
width reads 

r(j/v> - 00) = Pp - S 4 ( 8 7^ } E I a*5UI 2 (4-2) 

in the J/ip rest frame (M^ is the J/ip mass). If the J/ip is produced in e + e~ annihilations 
it is transversely polarized with respect to the beam direction and the angular distribution 
of the baryons emitted in the J/ip rest frame then exhibits a 1 + cos 2 9 dependence (up 
to corrections of (D(m Bs /M^)) characteristic of perturbative QCD fj]. The function p ps . 
in ( |4.2j) is the usual phase space factor 

Pp . s .(z) = VT^4^ • (4.3) 

As we said in the introduction we are going to calculate the invariant B Ba within the 
modified perturbative approach proposed in || , thus generalizing our analysis of the J/ i/j 
decay into nucleon-antinucleon pairs 0. As in previous perturbative calculations |TJ, ^ 
the J/ip meson is treated as a non-relativistic cc system and 0(v 2 /c 2 ) corrections are 



M 



Ai A2A 



11 2ni Bs 



neglected. In contrast to the decays of P-wave charmonium |L8| the ccg Fock state is 
suppressed in exclusive decays by inverse powers of the large scale which is provided by 
the c-quark mass m c [T^, pOfl , relative to the cc Fock state and is, therefore, neglected 
in our analysis. The use of m c , strictly speaking 2m c , as the large scale rather than the 
charmonium mass is consistent with the neglect of relativistic corrections. It is also well in 
the spirit of a perturbative approach since in the internal c-quark propagators the c-quark 
mass appears. The J/i/j state is written in a covariant fashion 



J/V; q, X) = M% )-=(i + M^(A) (4.4) 




where a and b are colour indices and is the J/ip decay constant being related to the 
J I ip wave function at the origin of the configuration space. Merely that part of the wave 
function is, to a reasonable approximation, required in the calculation of B since, as we 
already mentioned in the introduction, the c and the c quark only annihilate if their 
mutual distance is less than about l/m c [[| which is smaller than the J/ip radius [pi]] . 
is replaced by 2m c in the calculation of B%g and the baryon masses are neglected. Only 
the phase space factor in (|4.2|) is evaluated with the physical masses, and m^. 

The invariant B Bs receives its dominant contribution from the graphs with three in- 
termediate gluons, see Fig. [|. Within the modified perturbative approach the three-gluon 
contribution Bf g a to the J/ip decay into B 8 B 8 is of the form 

&ig = ^J[dx][dx']J ^^f H (x,x',b)ex V [-S(x,x',b,2m c )} (4.5) 
x (rr, b)*?i 3 (x', b) + \(^U^ b) + ^(x, b))(tf*(s', b) + *g\(x', b))~ 

for B$ = N, E, 5. This contribution is the same for all octet baryons belonging to the 
same isospin multiplet. Since the form of the A Fock state (2.1b) is slightly different the 
three-gluon contribution reads 

K = ^/[dx][da;']/^^f^(x,x',b)exp[-5(x,a;',b,2m c )] (4.6) 
x ^ 23 (x,b)^ 23 (x',b) + l(^ 23 (x,b) - ^ 21 (a;,b))(^ 2 3(x',b) - ^(x^b))' 

in the case of the A. The representation of B^ g as a convolution of wave functions and 
a hard scattering amplitude Th can formally be derived by using the methods described 
in detail by Botts and Sterman ||. The bj, canonically conjugated to the transverse 
momenta kj_j, are the quark separations in the transverse configuration space, bi and b 3 
correspond to the locations of quarks 1 and 3 in the transverse plane relative to quark 2 
and b 2 = b x — b 3 . ^| represents the Fourier transform of the wave function (see 
Sect. 2). 

Th is the Fourier transform of the usual momentum space hard scattering amplitude 
to be calculated from the Feynman graphs shown in Fig. |2]a. Up to corrections of order 





Figure 2: Feynman graphs for the J/ip decay ointo baryon-antibaryon. a) Three-gluon 
contribution (graphs with permutated gluon lines are not shown), b) electromagnetic 
contribution. 

oq, m 2 Bs /(4ml) and 5 2 /(4m 2 ) the hard scattering amplitude in b space reads 
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The quantities 



qi = 2[x l {l-x' i ) + {l-x l )x' l ] m 2 



~g\ = Ax^ m; 



(4.7) 



(4.t 



represent the virtualities of the internal quarks and gluons at zero transverse momenta. 
Since, as shown in |5[], the hard scattering amplitude only depends on the sum of transverse 
momenta, kj_j + k'jj, the transverse separation of any two quarks inside the baryon is the 
same as that of the corresponding antiquarks inside the antibaryon (b, = b'j). Physically, 
this property of the hard scattering amplitude means that baryon and antibaryon are 
created with identical transverse configurations of the quarks and antiquarks, respectively. 
The auxiliary variable b in ( [4-7]) serves as a Lagrange multiplier to the constraint Yl kj + 
k'j = 0. Since the virtualities of the gluons are timelike, T H includes complex-valued 
Hankel functions H that are related to the usual modified Bessel functions K , appearing 
for space-like propagators, by analytic continuation. 

The Sudakov factor exp[— S] entering ( |4.5| ) and ( [4.6| ) takes into account those gluonic 
radiative corrections not accounted for in the QCD evolution of the wave function as 
well as the renormalization group transformation from the factorization scale fip to the 
renormalization scales ti at which the hard amplitude Th is evaluated. The Sudakov 
factor, originally derived by Botts and Sterman || and later on slightly improved, can 
be found for instance in [22j| . The renormalization scales ti are defined in analogy to the 
case of electromagnetic form factors |§ as the maximum scale of either the longitudinal 
momentum or the inverse transverse separation associated with each of the gluons 



U = max(gi, g h l/b 3 ) , t 2 = max^, l/6 2 ) , ^3 = max(g 3 , g 3 , 



(4.9) 



Infrared cut-off parameters bi appear in the Sudakov factor which are naturally related 
to, but not uniquely determined by the mutual separations of the three quarks |23| . 




Figure 3: Percental accumulation of the three-gluon contribution to the width of the NN 
channel from regions of internal momenta where Yli=i a s(ti) < ( a s rit ) 3 - 

Following H] we chose hi = b = max{6i, b 2 , 63}. With this "MAX" prescription the three- 
gluon contribution B^ g 8 is unencumbered by as singularities in the soft end-point regions. 
As a consequence of the regularizing power of the "MAX" prescription, the perturbative 
contribution saturates in the sense that the results become insensitive to the inclusion of 
the soft regions. A saturation as strong as possible is a prerequisite for the self-consistency 
of the perturbative approach. The infrared cut-off b marks the interface betweeen the non- 
perturbative soft gluons, which are implicitly accounted for in the baryon wave function, 
and the contributions from soft gluons, incorporated in a perturbative way in the Sudakov 
factor. Obviously, the gliding factorization scale to be used in the evolution of the wave 
function, has to be chosen as \ip = 1/b. 

As an inspection of ( |4.5| ), (|4.6|) and ( |4.7| ) reveals, a nine dimensional numerical inte- 
gration has to be performed.0 Although this is a rather involved technical task it can be 
carried through with sufficient accuracy if some care is put into it. The numerical results 
are obtained from the wave functions discussed in Sect. 2 and for the following values 
of the J/ip decay constant and the c-quark mass: = 409 MeV, m c = 1.5 GeV. We 
evaluate a s in the one-loop approximation with rif = 4 and Aqcd — 210 MeV [21]. 



Before turning to the detailed discussion of the numerical results for the decay widths, 
we want to focus on an important feature of the modified perturbative approach. As 
mentioned in the introduction, one of the motivations for including the transverse hadronic 
structure and the Sudakov factor in the analysis is to achieve a theoretically self-consistent 
calculation in the sense that the bulk of the perturbative contribution is accumulated in 
regions where the strong coupling constant a s is sufficiently small. A method to check 
whether or not this is the case is to set the integrand in ( |4.5|) or ( [4.6|) equal to zero in 
those regions where n 3 =i a s{U) > ( a s nt ) 3 an d to evaluate the three-gluon contribution to 
the decay width as a function of af 11 . In Fig. |3| we show the accumulation profile for the 
nucleon case; it is typical of all baryons. Consulting Fig. || one sees that almost the entire 
result is accumulated in the comparatively narrow region of a s between 0.4 and 0.6. The 
regions with II 3 =i a s(ti) < 0.469 3 provide 50 % of the total result. Hence, our calculation 
of the J /ip decay widths into octet baryon-antibaryon pairs is theoretically self-consistent. 

2 Taking into account relativistic corrections to the J/ip wave function, i.e. its transverse momentum 
dependence, one would have to perform a 14 dimensional numerical integration which seems impossible 
with present day computers to a sufficient degree of accuracy. 





r 39 [eV] 


Data 


Channel 


Set 1 


Set 2 


Set 3 


Set 4 


0AS 


r exp [eV] 


PP 


174 


174 


174 


174 


140 


186 ± 14 




136 


128 


113 


108 


97.8 


110 ±15 


AA 


140 


133 


117 


107 


99.7 


117 ± 14 




107 


92.8 


62.5 


47.4 


46.9 


78 ± 18 



Table 3: Results for the decay widths of J /if) into octet baryon-antibaryon pairs for the 
four sets of distribution amplitudes defined by the expansion coefficients B n quoted in 
Tab. |2| (f^ = 409 MeV, m c = 1.5 GeV). For comparison we also quote the experimental 
results and, in the column labelled 0as, predictions evaluated with distribution amplitudes 
constructed from the asymptotic distribution amplitude instead from (|2.10|) (with m s = 
350 MeV). 



In Tab. |3] we present our results for the J /if) — > B$B S decay widths using the distri- 
bution amplitudes discussed in Sect. 2. For the sake of comparison we also expose the 
available experimental data |23jQ. As can be seen from the results obtained with the 



distribution amplitudes termed set 1 the phase space factor p p s is an important but not 
sufficient element for the suppression of the hyperon channels. Since m s = for these 
distribution amplitudes the differences in the predictions for the widths, except for the 
AA case, are only due to p p s .. For m s > the additional suppression of the end-point 
regions leads to smaller decay widths for the hyperon channels. In order to demonstrate 
the strength of that reduction we display the invariant i3^ 8 versus m s in Fig. |j. The 
number of strange quarks embodied in a given baryon is reflected in differently strong 
m s dependences of £>f^ 8 . As inspection of Tab. |3] brings to view, the phase space cor- 
rected three-gluon contributions nicely reproduce the experimental pattern of the decay 
widths provided they are computed with the distribution amplitudes of set 3. The value 
of 350 MeV for the parameter m s used in the construction of these distribution ampli- 
tudes appears reasonable, considering the interpretation of the mass factor ( |2.11j ) as the 
BHL exponential ||15| . Note that, as a consequence of the use of the c-quark mass in the 



calculation of Bf*, the result for the decay width of the pp channel given in Tab. |3] differs 
from that one reported in || . 

The distribution amplitudes defined by (|2.10| -2.12c) exhibit a little asymmetry with 
respect to permutations of the Xj. In order to study the importance of this asymmetry we 
also show in Tab. [3] results computed with the asymptotic distribution amplitude for the 
nucleon and with modifications of it for the hyperons that are constructed analogously 
to the distribution amplitudes described in Sect. 2. The asymmetry in the distribution 
amplitudes is seen to increase the magnitudes of the decay widths by about 30 % for the 
SS channel and about 20% for the other channels while the pattern of the predictions 
remains unchanged. 

Let us now assess the uncertainties of our predictions. The value of the J /if) decay 
constant used by us is determined from the leptonic J/tf> decay width. Since the J/if> 

3 We use PDG averages throughout. The original data are from pH EgL G% 
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Figure 4: The three-gluon contribution B^ g s vs. the mass parameter m s for the decay 



channels EE (solid), SS (dashed) and AA (dotted). 



decay constant, or more generally, the decay constant, f n , of a n 3 Si quarkonium state is 
defined by 



(0\jL\n 3 S l ) = f n M n 



the leptonic decay width of a n 3 Si state reads 



T(nS 



e e 



4?r e 2 Q a 2 fl 
3 M n 



(4.10) 



(4.11) 



where cq is the charge of the heavy quark the charmonium state consists of. Reexpressing 
( |4.11| ) in terms of the non-relativistic wave function at the origin of the configuration 
space, one arrives at the famous van Royen-Weisskopf width []29[| . The use of the decay 
constant instead of the wave function at the origin implies that we relate the J/tp — > BB 
(or the n 3 Si — > BB) widths to the leptonic width. By that means the uncertainties in the 
determination of the wave function at the origin via the usual van Royen-Weisskopf width 
cancel to a large extend. These uncertainties arise from relativistic and QCD corrections 
which seem to be large |30[ but are not well known 



21, 31 



The next uncertainty to be mentioned arises from the choice of the c-quark mass 
value. In accordance with calculations of the charmonium spectrum within non-relativistic 



potential approaches [^T[ and with a global fit of charmonium parameters |32| we take 
1.5 GeV as the favoured value. That value has, for instance, also been used in a recent 
analysis of P-wave charmonium decays into two pions [pl|. In spite of this, little changes 
of the m c value cannot be excluded and lead to an approximate rescaling of the decay 
widths by the factor (1.5 GeV/m c ) 8 . 

The value of Aqcd is also subject to uncertainties. A change of that value by, say, 
±20 MeV which roughly represents the inaccuracy of our present knowledge of Aqcd 



would alter the theoretical decay widths by about ±25%. We stress that for any changes 
of the m c and Aqcd values the ratios of any two decay widths calculated by us remain 
approximately unchanged. This assertion does not only refer to the J/ip — > B$B$ decay 
widths but it also applies to the still to be discussed J/i[) — > B 10 B W and ip' — > BB widths. 

In addition to the three-gluon contribution £>f^ 8 there is a subdominant, although in 
some cases perhaps sizeable, isospin- violating electromagnetic one, £>^ p3| , arising from 



the graph shown in Fig. |2|b. This contribution is proportional to the time-like magnetic 
form factor of the baryon. For the proton £>£ m amounts to about 15 % (in absolute value) 
of the total B p as is estimated from the recent measurement of the proton magnetic 
form factor in the time-like region [34]. Since the relative phase between B^ g and B p n 



is unknown we cannot simply add the two contributions 0. Therefore, we merely can 
state that, depending on the value of the relative phase, the electromagnetic contribution 
(including the interference term between it and B^ g ) to the J/tj) — > pp decay width can be 
as large as 30 % or only 2% . The comparison of our result with data indicates a relative 
phase close to ±7r/2 in the pp case. 

At this point a remark concerning the nn decay channel is in order. Since, as said 
repeatedly, the three-gluon contribution respects isospin symmetry any difference between 
the pp and nn decay widths must be due to the electromagnetic contribution. From 
experiment it is known that the widths for J/ip —>■ pp and J/ip —>■ nn decays agree within 
the experimental errors |^S| and that the time-like form factors of the proton and the 
neutron are approximately equal in modulus at least at s = 5.4 GeV 2 [j3~5]| . Thus, one 
may conclude that the relative phases between the three-gluon and the electromagnetic 
contributions are the same (up to a possible sign) for the proton and the neutron channel. 

The size of the electromagnetic contribution to the hyperon channels may be estimated 



from a recent analysis of the octet baryon form factors within a diquark model [36]. With 



the help of a few rather well determined parameters that model is able to describe a 
large number of exclusive observables. In particular relevant for the present work is the 
prediction that, in the space-like region, the magnetic form factors of the E + and E~ 
have opposite signs and are comparable in magnitude to that of the proton. The form 
factors of the A and E°, on the other hand, turn out to be very small. Predictions for 
the S~ form factor are not reported in [RBI but that form factor is presumably smaller in 
absolute value than the proton form factor. Assuming similar relative magnitudes of the 
form factors in the time-like region, we expect that, in the case of the hyperon channels 
listed in Tab. |3], the three-gluon contributions should match with the experimental data. 
This may not be the case for the E~E + channel; while the three-gluon contribution is the 
same as for the E°E° channel, the electromagnetic contribution may be large. 

There is also a small contribution to the invariant B Bs from cc annihilations mediated 
by two gluons and a photon. The gg'j contribution to B Bs being proportional to the ggg 
contribution [03 |, amounts to less than about 1% of the latter and is therefore neglected. 

Finally, from the measurement of the angular distribution of BgB 8 pairs produced in 
e + e~ — > J/ip — > B$B$ |26j one observes small violations of the helicity sum rule: The frac- 
tion of pp and AA pairs with equal helicities amounts to about 10% of the total number 
of pairs with, in particular in the A case, large errors. For the other hyperon channels the 
errors are so large that no conclusion can be drawn. The small amount of equal helicity 
pairs is what is to be expected if the process is dominated by perturbative QCD: Each 
of the virtual gluons creates a quark and an antiquark with opposite helicities. Since our 
baryon wave functions do not embody any non-zero orbital angular momentum compo- 
nent the quark helicities sum up to the baryon helicity. Hence, baryon and antibaryon 
are produced with opposite helicities. The small amount of B$B$ pairs with the wrong 
helicity combination observed experimentally, while indicating the presence of some soft 
contributions, can be considered as a hint that perturbative QCD is the dominant dy- 



The time-like nucleoli form factor is likely not under the regime of perturbative QCD in the energy 
region of interest; it is about a factor of 3 larger in absolute value than the form factor in the space-like 
f 2 



region at Q 2 = M3 



namical mechanism in the J/if> — > B$Bg decays. One should, however, be aware of these 
contributions when theoretical results for the J/if> — > B$Bg decays are compared with 
experiment. The production of B 8 B$ pairs with equal helicities can perhaps be explained 



as an constituent quark mass and/or baryon mass efffect [33. 37 



5 J/tp decays into decuplet baryons 

We are now going to apply the modified perturbative approach to the J /if) decays into 
the decuplet baryon-antibaryon channels A ++ A and £*"£*+ in the same manner as for 
the octet baryon case. As we said above, within a perturbative approach and with wave 
functions of zero orbital angular momentum components in the direction of the baryon 
momentum, B w and B w are produced with opposite helicities. For the same reason only 
helicities ±1/2 are possible. Thus, as for the octet baryon case, the only non-zero helicity 
amplitudes are the Ml™ x . They are fed by only one invariant B Bl ° (or, depending on the 
definitions of the covariants, only one linear combination of the invariants) out of the five 
the general covariant decomposition of the J /if} — ► B\oB w helicity amplitudes comprises. 
Within the modified perturbative approach, the three-gluon contribution to the invariant 
B Bl ° reads 




f d 2 b d 2 b- 

[dx][dx] J —^—-r^f H (x,x',b) exp[-S(x,x',2m c )} 



x *3g{x,b)*%Stf,b), (5.1) 

for Bio = A, £*. The three-gluon contribution is evaluated with the decuplet wave func- 
tions introduced in Sect. 3 and with the hard scattering amplitude ( [1.71) . Insertion of 
B B g° into ( |4.2j ) provides the wanted decay widths. It still remains to choose plausible 
values of the parameters a(io) and /(io) . The fact that we assume the same form for the 
k± dependence of the A and the nucleon wave functions, and that we use 4>as fc> r the A 
distribution amplitude which does not differ from the actual nucleon distribution ampli- 
tude fl2.10| ) greatly, suggests, as a first attempt, the ansatz a(io) = ct(8) (= 0.75 GeV -1 ) 



and P^ = P^. Thereby, it is perhaps plausible to evaluate the valence quark probability 
of the nucleon from 0as instead from ( |2.1(J| ). Doing so we find P^(AS) = 0.163 and the 
requirement P^ = P^(AS) leads to /(i ) = 0.0163 GeV 2 which is larger than /( 8 ) by a 
factor of \/6, i.e., the SU(6) result. The results for the J/ip decay widths into B w Bi 
pairs are presented in Tab. |j. As can be seen our result for the £*£* channel is too large 
as compared to the data while agreement is achieved for the AA case. 

Bearing in mind that the A is in completely symmetric flavour and spin states and that 
the Pauli principle effectively induces an additional repulsive interquark force, one may 
expect a larger radius, and hence a larger value of ano), for the A than for the nucleon. 
Therefore, we also try values for ano) slightly larger than that for ag 8 and fix in each 
case /(io) from the requirement P^t = 0.163 as before. In Tab. [| we list results obtained 
with the values 0.80 and 0.85 GeV -1 for a(io). For both these values of <2(io) satisfactory 
agreement with experiment is found. The uncertainties of the theoretical results are the 
same as for the octet baryons. In particular, one has to consider the possibility of large 
electromagnetic contributions (see Fig. §b). 



a (10) 

[GeV- 1 ] 


/(10) 

[GeV 2 ] 


Bi 


B 2 


B 3 


B 4 


B 5 


r£ A [eV] 


rgr [eV] 


0.75 


0.0163 


-0.494 


0.165 


-0.203 


-1.013 


0.058 


105 


66.1 


0.80 


0.0143 


-0.547 


0.182 


-0.216 


-1.081 


0.062 


82.6 


51.8 


0.85 


0.0127 


-0.601 


0.200 


-0.229 


-1.142 


0.065 


65.1 


40.8 




T exp [eV] 


96 ±26 


45 ±6 



Table 4: /(io), the expansion coefficients B n for the E* distribution amplitude and results 
for the three-gluon contribution to the J/ip — > A ++ A , decay widths for various 

values of the transverse size parameter O(io) (m s = 350 MeV). /no) is fixed by the 
requirement P^t = 0.163. 



6 Decays of the ip f and other quarkonia 

The extension of our approach to baryonic decays of the ip' (= ip(2 3 Si)) is now a simple 
matter. It is however important to realize that, in contrast to other authors |1[, we eval- 
uate the three-gluon contribution with the c-quark mass and not with the charmonium 
mass. This is, as we said, legitimate in a non-relativistic treatment of the charmonia. 
Hence, in order to get the ip' widths in our approach we have not to rescale the corre- 
sponding J/ip widths by (M^/M^/) 8 but rather by 

W - BB) = P^ m *l M j] ry^ e+ f) - BB) (6.1) 

1 ^ P P .sXm B /M^ r(J/^^e+e-) W ^ ; 1 ; 

which holds for both, octet and decuplet baryons. The charmonium decay constants are 
replaced by the leptonic decay widths by means of formula (|4.11|) . 

As an immediate examination of our approach one may apply the scaling relation ( |6.1| ) 
directly to the experimental data (see Tabs. j|, |], §). Considering the uncertainties due 
to the electromagnetic contributions, which in one or the other case may be large, ( |6.1| ) 
works quite well in particular for the octet baryons. For the decuplet baryons, on the other 
hand, it seems that the suppression of the ip' decay widths is a slightly underestimated, 
although the large experimental errors prevent any definite conclusion at present. One 
may suspect the neglect of the decuplet baryon masses to be responsible for that possible 
imperfection. In any case, an additional strong suppression, as provided by (M^/M^/) 8 
0, is in clear conflict with the data 0. This observation supports our attempt of using 
the c-quark mass in the calculation of the decay amplitudes rather than the mass of the 
charmonium state in question. 

Results for baryonic decay widths of the ip', evaluated through ( |6.1| ) from the set 3 J/ip 
widths, are listed in Tab. |5| where also recent experimental results of the BES collaboration 
are quoted. The data are still preliminary. The agreement between theoretical results 
(with a(io) — 0.85 GeV -1 in the decuplet baryon case) and experiment is generally good 
although our results seem to be a bit too large for the EE and channels. For a 

discussion of uncertainties we refer to Sect. 3. 

Computation of the ^(S 3 ^) decay widths are difficult within our approach. The 
relativistic corrections are presumably larger since the ^(3 3 Si) mass is above the threshold 
for open charm production. The BB decay widths are likely to be very small and it is 



channel 


pp 




AA 




A++A— 






76.8 


55.0 


54.6 


33.9 


32.1 


24.4 


Texp I 

m 


76 ± 14 
53 ± 15 


26 ± 14 


58 ± 12 


23 ±9 


25 ±8 


16 ±8 



Table 5: The three-gluon contribution to the if)' — > 55 decay widths (in [eV]) computed 
through (|6.1|) from the set 3 J/ip widths (see Tab. ||). a(io) = 0.85 GeV -1 . 



hardly conceivable that they will be measured. Therefore, we refrain from estimating 
these decay widths. 

Results for bottonium decays, on the other hand, can safely be calculated within our 
approach. The hard scale, provided by the 6-quark mass, is larger than in the charmonium 
case and relativistic corrections are smaller. But, as it turns out, the predicted decay 
widths for the baryonic channels are also very small. Approximately, i.e. ignoring the fact 
that the fcj_-dependent suppression of the three-gluon contribution is a little bit different 
in the two cases, we find the following rescaling formula 



r(T - BB) = 4 PpA r B '^\ T 5 ( T Tr e+ !"\ (- ) n W - BB) 
1 ' P P .s.(m B /M^ r(J/^e+e-) \mj 



(6.2) 



Using m b = 4.5 GeV we obtain, for instance, a value of 0.03 eV for the T — > pp decay width 
which value corresponds to a branching ratio of 0.5 x 10~ 6 well below the experimental 
upper bound p4|. The decay widths for the other BB channels are even smaller. 



7 Summary 

In this investigation, we applied the modified perturbative approach to the decays of J j ip 
and ip' into baryon-antibaryon pairs. We demonstrated that, on the basis of plausible 
baryon wave functions for which SU(3)f symmetry is only mildly broken by quark mass 
effects and for which even SU(6) symmetry, in sharp contrast to the QCD sum rule based 
wave functions, approximately holds, the experimental data for octet and decuplet baryon 
channels are quite well reproduced by the phase space corrected three-gluon contributions. 
The perturbative contributions to the decay widths are calculated self-consistently in the 
sense that the bulk of a perturbative contribution is accumulated in regions of reasonably 
small values of a s . Besides the form of the wave functions used by us our analysis differs 
from previous ones in the following points: 

i) The use of the modified perturbative approach allows to take into account the running 
a s and the evolution of the wave functions properly, in contrast to the usual leading 
twist analysis. The virtualities of the internal c quarks and gluons can be chosen as the 
arguments of a s . Hence, a s reflects the characteristic scale of the process under question. 
The running coupling constant is therefore not a quasi-free parameter that can, within 
a certain range, be chosen arbitrarily. Since the decay widths are proportional to a^, a 
large factor of uncertainty is therefore hidden in the standard perturbative analysis. 

ii) The hard scattering amplitude is computed with the c-quark mass instead of the 



charmonium mass. The use of the c-quark mass is consistent with the non-relativistic 
treatment of the charmonium state and with the perturbative approach. Hence the ip' - 
J I ip scaling relation ( |6.1| ) holds in our approach. That scaling relation is nicely confirmed 
by the data in contrast to the usual (M^/M^) 8 scaling. 

Although our results for the decay widths agree with the data, we are aware of a 
number of uncertainties in our calculation. Above all, we mention as a source of uncer- 
tainty the value of the c-quark mass. The value we have chosen (1.5 GeV) is consistent 
with other constraints on that mass. Another uncertainty arises from the electromagnetic 
contribution which, at least in the pp case, can be large. In virtue of the unknown relative 
phase between the electromagnetic and the three-gluon contribution, the first cannot be 
taken into account properly. We emphasize that most of the uncertainties cancel in ratios 
of widths to a large extend. 

It should be noted that we have not found hints at substantial contributions from the 
ccg Fock state, i.e., colour octet contributions fI8 |, neither in the baryonic J/ip nor in the 
ip' decays. 

An interesting class of J/ip decays are the B^B\q channels. While the three-gluon 
contributions to the pA and E*°A channels are strictly zero |38| they are non-zero - 
although small - for other BgB w channels, since our wave functions exhibit only a mild 
breaking of SU(3)f symmetry. Experimentally, only upper bounds are known for the first 
two channels |24| saying that these decays are indeed suppressed by at least an order of 



magnitude as compared to, say, the pp or the AA channels. The experimental widths 



for the other B 8 Biq channels are surprisingly large | 24]| . Thus, in accordance with |38|j , 
we expect as the dominant SU(3)f breaking mechanism for these reactions a sizeable 
electromagnetic contribution. 

Finally, we have to mention that there are a few exclusive charmonium decays which 
cannot be described within the standard or the modified perturbative approach. Thus, 
the relatively large branching ratio of the process J/ip — > pit observed experimentally, 
for instance, indicates a substantial violation of hadronic helicity conservation while only 
mild violations are observed in the baryonic J/ip decays (see the discussion in Sect. 3). 
For a discussion of this puzzle and a possible solution of it by means of a hypothetical 



glueball see |39 . 
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